Search for: All records

Models currently fail to reproduce observations of the coldest regions in the Sun’s atmosphere, though recent work suggests the thermal Farley–Buneman instability (TFBI) may play a critical role. This meter-scale, electrostatic, multifluid plasma instability causes turbulence and heating in the coldest regions of the solar chromosphere. This paper describes how TFBI turbulence and heating varies across multifluid 2D, kinetic 2D, and kinetic 3D simulations. It also presents the first 3D simulations of the TFBI. We find that multifluid and kinetic 2D simulations produce similar results overall, despite using vastly different approaches. Additionally, our kinetic 3D simulations produce a similar or somewhat larger amount of heating compared to 2D, as contributions from the parallel electric field account for only (13 ± 2.5)% of the total turbulent heating in 3D. 

The ponderomotive force has been suggested to be the main mechanism to produce the so-called first ionization potential (FIP) effect—the enrichment of low-FIP elements observed in the outer solar atmosphere, in the solar wind, and in solar energetic events. It is well known that the ionization of these elements occurs within the chromosphere. Therefore, this phenomenon is intimately tied to the plasma state in the chromosphere and the corona. For this study, we combine IRIS observations, a single-fluid 2.5D radiative magnetohydrodynamics (MHD) model of the solar atmosphere, including ion–neutral interaction effects and nonequilibrium (NEQ) ionization effects, and a novel multifluid multispecies numerical model (based on the Ebysus code). Nonthermal velocities of Si IV measured from IRIS spectra can provide an upper limit for the strength of any high-frequency Alfvén waves. With the single-fluid model, we investigate the possible impact of NEQ ionization within the region where the FIP may occur, as well as the plasma properties in those regions. These models suggest that regions with strongly enhanced network and type II spicules are possible sites of large ponderomotive forces. We use the plasma properties of the single-fluid MHD model and the IRIS observations to initialize our multifluid models to investigate the multifluid effects on the ponderomotive force associated with Alfvén waves. Our multifluid analysis reveals that collisions and NEQ ionization effects dramatically impact the behavior of the ponderomotive force in the chromosphere, and existing theories may need to be revisited. 

Abstract Models fail to reproduce observations of the coldest parts of the Sun’s atmosphere, where interactions between multiple ionized and neutral species prevent an accurate MHD representation. This paper argues that a meter-scale electrostatic plasma instability develops in these regions and causes heating. We refer to this instability as the Thermal Farley–Buneman Instability (TFBI). Using parameters from a 2.5D radiative MHD Bifrost simulation, we show that the TFBI develops in many of the colder regions in the chromosphere. This paper also presents the first multifluid simulation of the TFBI and validates this new result by demonstrating close agreement with theory during the linear regime. The simulation eventually develops turbulence, and we characterize the resulting wave-driven heating, plasma transport, and turbulent motions. These results all contend that the effects of the TFBI contribute to the discrepancies between solar observations and radiative MHD models. 

Context. Spectroscopic observations of the emission lines formed in the solar transition region commonly show persistent downflows on the order of 10−15 km s −1 . The cause of such downflows, however, is still not fully clear and has remained a matter of debate. Aims. We aim to understand the cause of such downflows by studying the coronal and transition region responses to the recently reported chromospheric downflowing rapid redshifted excursions (RREs) and their impact on the heating of the solar atmosphere. Methods. We have used two sets of coordinated data from the Swedish 1 m Solar Telescope, the Interface Region Imaging Spectrograph, and the Solar Dynamics Observatory for analyzing the response of the downflowing RREs in the transition region and corona. To provide theoretical support, we use an already existing 2.5D magnetohydrodynamic simulation of spicules performed with the Bifrost code. Results. We find ample occurrences of downflowing RREs and show several examples of their spatio-temporal evolution, sampling multiple wavelength channels ranging from the cooler chromospheric to the hotter coronal channels. These downflowing features are thought to be likely associated with the returning components of the previously heated spicular plasma. Furthermore, the transition region Doppler shifts associated with them are close to the average redshifts observed in this region, which further implies that these flows could (partly) be responsible for the persistent downflows observed in the transition region. We also propose two mechanisms – (i) a typical upflow followed by a downflow and (ii) downflows along a loop –from the perspective of a numerical simulation that could explain the ubiquitous occurrence of such downflows. A detailed comparison between the synthetic and observed spectral characteristics reveals a distinctive match and further suggests an impact on the heating of the solar atmosphere. Conclusions. We present evidence that suggests that at least some of the downflowing RREs are the chromospheric counterparts of the transition region and lower coronal downflows.